Plasma wake inhibition at the collision of two laser pulses in an underdense plasma 
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An electron injector concept for laser-plasma accelerator was developed in ref pj and Q] ; it 
relies on the use of counter-propagating ultrashort laser pulses. In Q], the scheme is as follows: 
the pump laser pulse generates a large amplitude laser wakefield (plasma wave). The counter- 
propagating injection pulse interferes with the pump laser pulse to generate a heatwave pattern. 
The ponderomotive force of the heatwave is able to inject plasma electrons into the wakefield. We 
have studied this injection scheme using ID Particle in Cell (PIC) simulations. The simulations 
reveal phenomena and important physical processes that were not taken into account in previous 
models. In particular, at the collision of the laser pulses, most plasma electrons are trapped in the 
heatwave pattern and cannot contribute to the collective oscillation supporting the plasma wave. At 
this point, the fluid approximation fails and the plasma wake is strongly inhibited. Consequently, 
the injected charge is reduced by one order of magnitude compared to the predictions from previous 
models. 



C*2 

o ■ 

^ ■ 



> 
o 

o 
m 
o 
t> 
o 

o 

>>: 

Oh- 



X 



An intense laser pulse can drive an electrostatic plasma 
wave via the ponderomotive force which scales as Va 2 , 
where a is the normalized potential vector of the laser 
: |a| = a = 8.6 x 10- 10 A[^m]/ 1/2 [W.cm- 2 ] for a lin- 
early polarized laser. When the pulse duration is close 
to the plasma period (A p /c), the laser pulse ponderomo- 
tive force pushes electrons and efficiently creates charge 
separation (ions hardly move). It results in a travelling 
longitudinal wave whose phase velocity v p is equal to the 
group velocity of the laser. In an underdense plasma, v p 
is very close to c, the speed of light, thus enabling accel- 
eration of electrons to very high energies once they are 
trapped in the wake [3[. But in a linear or moderately 
nonlinear regime, an electron with no initial velocity is 
not trapped by this travelling wave and consequently, not 
accelerated. In a more nonlinear regime, transverse wave 
breaking effects can result in the self-trapping of electrons 
in the so-called "bubble regime" J4]. This phenomenon 
has been observed in 2004 in @, 0] where quasi mo- 
noenergetic electron beams at the 100 MeV level were 
obtained. Nevertheless, in this scheme, self-injection and 
acceleration depend on the precise evolution of the laser 
pulse. Therefore, a fine control over the output electron 
beam is hard to achieve. On the contrary, precise control 
of electron injection would translate into good tailoring of 
the electron beam parameters, and would be most useful 
for applications (aOJ. 

To trap electrons in a regime where self-trapping does 
not occur, one has to externally inject electrons in the 
plasma wave, i.e. give electrons an initial momentum. 
In addition, electrons should be injected in a short time 
(< Ap/c) in order to produce a monoenergetic beam. 
This can be achieved using additional ultrashort laser 
pulses whose purpose is only restricted to triggering elec- 
tron injection. Umstadter et al. [l(| first proposed to use 
a second laser pulse propagating perpendicular to the 
pump laser pulse. The idea was to use the radial pon- 



deromotive kick of the second pulse to inject electrons. 
Esarey et al. [l[ proposed a counter-propagating geome- 
try based on the use of three laser pulses. This idea was 
further developed in Ref. 0] , where only two laser pulses 
are necessary. In this scheme, a main pulse (pump pulse) 
with maximum amplitude ao creates a high amplitude 
plasma wave and collides with a secondary pulse of lower 
maximum amplitude a\. The interference of the two 
beams creates a ponderomotive beatwave pattern with 
phase velocity v^w = 0, and thus enables to preaccel- 
erate background electrons. Upon interacting with this 
field pattern, some background electrons gain enough mo- 
mentum to be trapped in the main plasma wave and then 
accelerated to high energies. The force associated with 
this ponderomotive beatwave scales as Fb w — 2koaoo,\, 
where kg is the central wavenumber of both pulses, Fb w 
is many times greater than the ponderomotive force as- 
sociated with the pump laser -Fpond ~ k p a^ since in an 
underdense plasma fco >> k p . Therefore, the mechanism 
is still efficient even for modest values of ao and a\. As 
the overlapping of the lasers is short in time, the elec- 
trons are injected in a very short distance and can be 
accelerated to an almost monoenergetic beam. This con- 
cept has been recently validated in an experiment pH . 
using two counter-propagating pulses. Each pulse had a 
duration of 30 fs at full width half maximum (FWHM) , 
with ao = 1.3, a\ — 0.4. They were propagated in a 
plasma with electron density n e = 7 x 10 18 cto~ 3 corre- 
sponding to 7 P = ko/k p = 15. It was shown that the 
collision of the two lasers could lead to the generation of 
stable quasi-monoenergetic electron beams. The beam 
energy could be tuned by changing the collision position 
in the plasma. 

The precise understanding of these experiments, as 
well as the optimization of this process, motivate the 
present study. We have used ID Particle in cell (PIC) 
simulations to model electron injection in the plasma 
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wave at the collision of the two lasers, and their sub- 
sequent acceleration. The PIC simulations are compared 
to existing models [l[ and show major differences, such 
as the plasma fields behavior and the amount of injected 
charge. 

We first describe the fluid model developped in Ref. 
0, • In the linear approximation the wakefield poten- 
tial (due to charge separation) is a superposition of three 
potentials [l]| : $ = <f> + $i + &b- $o,i are the charge 
separation potentials driven by the laser pulses a 01 , $ 
representing the main accelerating structure we want to 
inject electrons in, and is driven by the beatwave. The 
normalized expressions of the wakefields fa = e$ /mc are 
given by : 
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where £o,i — z T v p t stands for the phase relative to 
the pump and injection lasers, k p is the plasma wavevec- 
tor and brackets < . > denote the time-average over 
the fast varying scale (l/(fcoc)). The term < ao.ai > 
is the beatwave, appearing only during the collision of 
the two lasers. Its spatial scale is given by 2kg, for ex- 
ample for circularly polarized lasers we have ao.ai = 
ao(r, £,o)ai(r, £i) cos(2k z). 

The first approximation of this analytical model con- 
sists of neglecting the last term, fa. This is supported 
by the fact that the density variation linked with this 
electrostatic potential scales as Srib/n — Ak^/k 2 fa with 
ko/kp » 1. This density variation being limited to 
5rib/n ^2 1 even in a non- linear regime, the potential can 
not be greater than fa ^2 k 2 /4/cq which is often negligible 
compared to the other terms scaling as 0o,i — &o l- 

The second assumption is that we can separate the 
dynamics of the particles inside and outside the beat- 
wave because the timescales are different (i.e. an elec- 
tron will see a constant wakefield during an oscillation 
in a beatwave bucket). An underlying hypothesis here 
is that there are two different species of electrons, those 
maintaining the wakefield, or fluid electrons, and those 
being trapped in the beatwave, or test electrons. 

Starting from these hypothesis, one can even build an 
analytical model when the lasers have modest intensities 
(ao < 1) and are circularly polarized. In that case, the 
beatwave pattern is time independent. Using this an- 
alytical model, one can find an analytical threshold for 
trapping ([Ell]) : 2^a ai > u z (£ mm) where u z (£_ Qmin ) is 
the minimum normalized longitudinal momentum p z /mc 
for which the electrons can be t rapp ed in the wakefield. 
The experiments described in ref[ll| operated well above 
this threshold. 

A numerical implementation of this model is a par- 
ticle tracking code @ where test particles are pushed 
in prescribed fields. The plasma fields are given by 
E p i = — V($o + $1) and the laser fields are solutions 



of the paraxial wave equation with a linear plasma re- 
sponse (gaussian beams). As expected, this model shows 
good agreement with the analytical model in the circu- 
lar polarization case and allows to extend the scheme to 
the linear polarization case for which analytical theory 
is untractable. For linear polarization, the beatwave has 
a fast var ying time dependence, which leads to stochas- 
tic effects [13}. In that case, the trapping thresholds are 
even lower. 

The main results of this model, that we will later on re- 
fer to as the "prescribed fields" model, are that electrons 
can be injected with modest values of ao and a\ (linear 
regime) . The resulting bunch is quasi- monoenergetic be- 
cause all the electrons are injected in a short distance 
(where the beating occurs). Moreover, the charge is ex- 
pected to be up to some hundreds of pC in the linear 
regime. ID PIC simulations, where fields are linked to 
the motion of electrons, have also already been carried 
out [14j . They confirmed that an electron beam with low 
energy spread, low emittance and short bunch length can 
be obtained with modest values of ao and a\. 

Here we compare ID PIC simulations with a ID pre- 
scribed field model. In this later model, we have used 
the ID nonlinear theory of wakefield generation to be 
able to compare the results in a more nonlinear regime 
(1 < ao < 2). The corresponding equation writes [la ]: 
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(b) t = -13 fs 
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FIG. 1: Longitudinal electric field computed at different times 
in PIC simulation (thick solid line), and in prescribed fields 
simulation (dotted line). 

ID PIC simulations were carried out with the code 
CALDER [l6(. For all the simulation results presented 
here, we have used a simulation box of 6400 cells, each cell 
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measuring 0.1/fco- Both lasers were gaussian pulses with 
30/s duration at FWHM and wavelength Ao = 0.8/im. 
To avoid edge effects, the main pulse entered the plasma 
through a density ramp of 100/im, the maximum density 
being 7 x 10 18 cm~ 3 . The simulation box was kept fixed 
to let the laser pulses enter the simulation correctly and 
we otherwise used a moving window to follow the main 
pulse over long propagation distances.. 

Fig. [T] shows the time evolution of the longitudinal 
electric field in ID PIC simulations when the polariza- 
tions of the two pulses are parallel (solid line), compared 
to the fields externally specified by equation (dot- 
ted line), for ao = 2 and a\ = 0.4. The electric field is 
given in linear wave breaking field unit E p — cm e ujpj ' e 
corresponding here to E p = 250GV/m. The pump pulse 
propagates from left to right. We have also plotted the 
transverse electric field (thin dotted line) to show the po- 
sition of the laser pulses. The collision occurs at z = 
at time t = 0. Fig[TJa shows the electric field 40 fs before 
the collision, it fits very well with the solution of equa- 
tion ([3|). During the collision of the two beams (Fig[T]b) 
we clearly see a small spatial scale pattern at 2kg, cre- 
ated by the beatwave. However, in this non-linear regime 
(Srib/n n 1) we do not observe a superposition of this 
pattern with the usual wakefields as is the case in ref 
The most remarkable feature is the strong distor- 
sion of the plasma wave at the position where the two 
pulses collide. This distorsion remains after a long time 
(FigfTJc, [TJd) and a numerical estimate shows that the 
plasma wave amplitude is decreased by a factor of 10 at 
this position. 
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FIG. 2: Typical trajectories of electrons in phase space, a): 
Pump laser only, linear polarization (ao = 2); b): Pump and 
injection lasers, parallel linear polarizations (ao = 2 and ai = 
0.4). c): Pump laser only, circular polarization (ao = 2/\/2); 
d): Pump and injection lasers, circular polarizations (ao = 
2/S/2 and ai = 0.4/\/2) 

This decrease of the wakefield amplitude follows from 
the fact that the electrons experiencing the collision of 
the lasers are trapped in the beatwave buckets and hence 
do not respond coherently (fluid-like) to the ponderomo- 



tive force of the main pulse. Without this fluid hypoth- 
esis, equation ([3]) fails. It would still hold if only some 
electrons were trapped, but here, as the phase velocity of 
the beatwave is zero, all electrons are trapped and kinetic 
effects cannot be treated as a perturbation. 

To better understand this plasma wave " damping" , we 
have performed simulations where test electrons are sub- 
mitted to the laser fields in vacuum. In particular, this 
allows to understand electron motion in the pondcromo- 
tive beatwave. Results of these simulations are plotted 
in figure [5] where we represent the orbits of test electrons 
in (z — ZQ,p z ) phase space. Electrons are initially taken 
at rest, their initial position zq is randomly distributed 
around the collision position. Fig. [5]a illustrates the tra- 
jectories of electrons experiencing the pump laser field 
only (with gq = 2) : in the laser reference frame, every 
electron has the same motion because they are pushed by 
the same ponderomotive force. This fluid motion leads 
to the charge separation which drives the electrostatic 
field known as the plasma wave. On the contrary, when 
a parallel polarized injection beam is included, trajec- 
tories are drastically modified (FigJ^lb). Electrons are 
locally trapped in beatwave buckets, there is no large 
scale charge separation and hence the plasma wake is no 
longer excited. This wake inhibition phenomenon also 
occurs for circularly polarized lasers as shown in FigsHJc 
andOd. This demonstrates that the plasma wave distor- 
sion is not related to the stochastic motion of electrons, 
but really to their trapping in the beatwave buckets. 

We now investigate the consequences of this plasma 
wake inhibition on the features of the electron beam. Fig- 
ure [3] shows the energy of the electron beam versus their 
phase with respect to the laser (main plots) as well as the 
corresponding spectra (inserts) 300 /mi after the collision 
of the two beams. FigJ3]a corresponds to the prescribed 
field simulation and Fig[3]b corresponds to the ID PIC 
simulation. Qualitative differences on the energy spectra 
are only minor, both spectra showing a peaked distri- 
bution around 45 MeV. Electrons are trapped mainly in 
the first wakefield bucket, and less than ten percent are 
trapped in the next ones. However, depending on the 
model, the trapped charge differs by almost one order of 
magnitude. 

In the prescribed fields model, electrons are preaccel- 
erated in the beatwave. As their initial velocity is lower 
than the wakefield phase velocity, they slip backward in 
phase. They are then trapped in the wakefield provided 
the energy gained in the beatwave was sufficient. In that 
case, electrons catch up with the plasma wave and are 
accelerated to high energies. In the more realistic PIC 
simulation, the beatwave preaccelerates electrons in the 
same way but the wakefield is inhibited and most elec- 
trons just slip back in phase space without being trapped. 
In order to be trapped, electrons need to be preacceler- 
ated at the end of the collision. As they slip back slower 
in phase space than fluid electrons, they witness a re- 
stored wakefield suitable for trapping and acceleration 
when they reach the back of the pump pulse. Therefore, 
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FIG. 3: Energy gain after 300/im as a function of relative 
phase £ = z — ct, (£ = corresponding to the maximum of 
the pump pulse). Top : prescribed fields simulation, bottom 
: PIC simulation. 



the plasma wake inhibition reduces the phase space vol- 
ume of injected particles. This leads to lower trapped 
charge and to a lesser extent, to a smaller energy spread 
in the PIC simulation. 

Figure [4] shows a more complete comparison of the 
charge obtained in both simulations (logarithmic scale), 
with ao varying beetween ao = 0.9 and ao = 2. Here, 
ao = 0.9 corresponds to the trapping threshold. The 
hollow squares represent the charge obtained in ID PIC 
simulations ; the triangles represent the charge obtained 
using the prescribed fields simulations. As we can see, the 



prescribed fields simulation overestimates the charge by 
a constant factor 8. The plasma wake inhibition and the 
reduction of the trapped charge has also been confirmed 
in 3D hybrid PIC simulations. [l7l ]. 
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FIG. 4: Comparison of injected charge obtained in PIC sim- 
ulations and in prescribed fields simulations. 

In conclusion, although the previous modelling of the 
collision of two laser pulses in an underdense plasma gives 
a good understanding of the general processes at play in 
the injection of electrons, we have observed and explained 
an important physical process neglected to date. Using 
ID PIC simulations, we have shown that the beatwave 
does not only preaccelerate electrons but also reduces the 
wakefield amplitude. In this process, the trapped charge 
decreases by an order of magnitude and to a lesser extent, 
the energy spread can be improved. Further studies will 
focus on minimizing this effect and increasing the injected 
charge. 
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